QCD@Work 2003 - International Workshop on QCD, Conversano, Italy, 14-18 June 2003 



Potential models, sum rules, color singlet model versus NREFT in Heavy 
Quarkonium 
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I briefly review how nonrelativistic effective field theories solve old puzzles and open problems in heavy quarkonium physics. 



1 Historical recollection 

Today quarkonium is a key issue in most of the accelera- 
tor experiments (see [[DEI) where millions and in some 
cases billions of heavy quarkonium states are being pro- 
duced. Historically quarkonium has been one of the most 
important playgrounds for our understanding of the strong 
interactions. 

In perspective, the seventies were dominated by the J/if/ 
discovery, around 3095 MeV with a lifetime about 1000 
times longer than that of other particles of comparable 
mass. This was known in the particle physics community 
as the November revolution. During the one year after the 
discovery more than seven hundreds papers were written 
related to Jj^i. Most of the papers dealt with quark and 
antiquark in a colour singlet state bound by a phenomeno- 
logical potential. Since the structure of the energy levels of 
charmonium and the soonly after discovered bottomonium 
is somehow intermediate bewteen a Coulomb and an har- 
monic oscillator structure and due to the asymptotic free- 
dom idea [ 3 1, the basic form of the static potential is taken 
as a superposition of a Coulomb and a linear potential: 
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where a s , the strong coupling constant at some scale, and 
<x, the string tension, have to be fit on the data. Such 
static singlet potential, together with similar phenomeno- 
logical spin-dependent corrections, when used inside a 
Schrodinger equation gives an overall successfull descrip- 
tion of heavy quarkonium spectra and decays [ 0] E] 0. 
However, already in these pioneering papers the problem 
arised of the relation of these potential models to QCD and 
to the QCD parameters. Starting with the work of Wilson, 
an effort [ 6 7 1 was devoted to relate the singlet potential, 
static and relativistic corrections, to QCD average values 
of Wilson loop and field strength insertions into the Wil- 
son loop, which are objects easy to evaluate in lattice QCD 
[|SJ. However, it turned out that these results for the poten- 
tials in terms of Wilson loops were inconsistent with the 
one loop calculation of the potentials in perturbation the- 
ory [Q. 



The color singlet potential model encountered soon con- 
crete problems also in the calculation of inclusive annihi- 
lation decays rates of heavy quarkonium states into light 
hadrons (hadronic) and photons and lepton pairs (elec- 
tromagnetic). It was assumed that the decay rate of the 
quarkonium state factored into a short distance part /, cal- 
culated in perturbative QCD as the annihilation rate of the 
heavy quark and antiquark, and a long distance nonpertur- 
bative part given in terms of the quarkonium wave function 
(or its derivatives) evaluated at the origin: 



T = /(a s (m)) • |.A(0)| 2 



(2) 



Explicit calculations at next to leading order in a s in per- 
turbation theory for S and for P wave decays into photons 
supported the factorization assumption which could not, 
however, be proved on general grounds for higher orders of 
perturbation theory. Indeed, in the case of P wave decays 
into light hadrons, it turned out that at order the fac- 
torization was spoiled by logarithmic infrared divergences. 
The same problem appeared in relativistic corrections to 
the annihilation decays of S wave states [|9]- Therefore, 
potential models were unable to supply at higher order an 
infrared finite prediction for the inclusive decays. 

The eighties witnessed the success of the sum rules. It was 
a distinctive prediction of the sum rules that as a result of 
QCD the rj c mass has to be located around 3 GeV and not 
at the much smaller value of 2.83 GeV claimed at that time 
by the experiments f II 01 - The later discovery of the gen- 
uine r\ c state with mass 2.98 Gev at Stanford was a great 
success of QCD and gave motivations for further work on 
the sum rules. Sum rules work in terms of Wilson opera- 
tor expansion. The Green function of the quark-antiquark 
pair injected in the vacuum is expanded in terms of per- 
turbative coefficients and local nonperturbative objects like 
the gluon condensate and the quark condensate. Physi- 
cally as far as the QQ distance r is held smaller than the 
confinement scale A.q CD the binding bewteen the quark is 
perturbative and the propagation of the heavy pair is tak- 
ing place in nonperturbative "external" vacuum fields. In 
such situation the interaction with the nonperturbative (at 
the scale Aqcd) vacuum gluonic fields can be expanded 



2 



QCD@Work 2003 - International Workshop on QCD, Conversano, Italy, 14-18 June 2003 



in multipoles and the leading contribution of the vacuum 
fields to the quarkonium energy levels is proportional to 
the gluon condensate <0|F£ v (0)F° y (0)|0) [112. Such con- 
densate corrections are incompatible with the genuine po- 
tential described in the previous paragraph. Precisely, the 
effects of the nonperturbative fluctuations of the gluonic 
field cannot be expressed in terms of a singlet QQ interac- 
tion potential and one should consider states that contain 
both the QQ pair and a gluonic excitation, including the 
case in which both of them are in a color octet \ II II . In 
brief, sum rules results for heavy quarkonia seemed to be 
in strong contradiction with the potential model picture. 

The nineties were dominated by the data on prompt pro- 
duction of charmonium at Fermilab. The first hadron col- 
lider measurements of inclusive charmonium production at 
CERN and from the CDF collaboration at Fermilab could 
not separate charmonium produced in hard scattering re- 
actions from charmonium produced in weak decays of B 
mesons. Thus the comparison with theory was uncer- 
tain. A rigorous test of the colour singlet production model 
(i.e. the assumption that the two quark were produced in 
a color singlet state) became possible with the CDF data 
on direct charmonium production where the contributions 
from the B-decays had been removed using microvertex- 
detection. With these data it became clear that the colour- 
singlet model failed dramatically when confronted with the 
experimental results r lT2"l . 

This historical excursion ends with several puzzles and 
open problems related to the existence or nonexistence of 
a quark antiquark singlet potential, and the validity or not 
of a color singlet picture. Today the stage is taken by the 
many experiments accumulating high statistic data sam- 
ples at quarkonium resonances and with several produc- 
tion mechanisms. It becomes therefore even more relevant 
to clarify the open theoretical problems and to supply a 
clean and under control QCD picture of these systems. In 
the next section we will show how nonrelativistic effective 
field theories provide the solution. 

2 Non Relativistic Effective Field Theories 
for Heavy Quarkonium 

The reason for which the EFT approach is so successfull 
for heavy quarkonium is the fact that heavy quarkonium, 
being a non-relativistic bound state, is characterized by 
a hierarchy of energy scales m, mv and mv 2 , with m the 
heavy-quark mass and v <K 1 the relative heavy-quark ve- 
locity. A hierarchy of EFTs may be constructed by sys- 
tematically integrating out modes associated to the energy 
scales not relevant for the quarkonium system. Such inte- 
gration is made in a matching procedure that enforces the 
complete equivalence between QCD and the EFT at a given 
order of the expansion in v and a s . 

Integrating out degrees of freedom of energy m, which 



for heavy quarks can be done perturbatively, leads to non- 
relativistic QCD rNROCDirfT3lfT4l. This EFT still con- 
tains the lower energy scales as dynamical degrees of free- 
dom. In the last years, the problem of integrating out the re- 
maining dynamical scales of NRQCD has been addressed 
by several groups and has reached a good level of under- 
standing (a list of references can be found in [ 18 1). The 
EFT obtained by subsequent matchings from QCD, where 
only the lightest degrees of freedom of energy mv 2 are left 
dynamical, is potential NRQCD, pNRQCD [QUEl- ™ s 
EFT is close to a quantum-mechanical description of the 
bound system and, therefore, as simple. It has been sys- 
tematically explored in the dynamical regime Aqcd < fnv 2 
in [ [n] [19 20 1 and in the regime mv 2 <k Aqcd < mv in [ 
I17ll21ll22l . The quantity Aqcd stands for the generic scale 
of non-perturbative physics. 

Inside the EFT, the power counting in the small quantity v 
enables to select the operators that contribute to physical 
quantities up to a definite order in v. The EFT approach 
makes it possible, in the case of several observables, to 
achieve a rigorous factorization between the high-energy 
dynamics encoded into matching coefficients calculable in 
perturbation theory and the non-perturbative QCD dynam- 
ics encoded into few well-defined nonperturbative contri- 
butions to be fitted on the data or calculated on the lattice. 
Thus, several model independent QCD predictions become 
possible. I will detail these in the following section for 
NRQCD and pNRQCD. 

3 Non Relativistic QCD 

NRQCD is the EFT obtained by integrating out the hard 
scale m. The mass m being larger than the scale of 
non-perturbative physics, Aqcd, the matching to NRQCD 
can be done order by order in a s . Hence, the NRQCD 
Lagrangian can be written as a sum of terms like 
fn Olf" '/ 'm d, <~ A ', ordered in powers of a s and v. More specif- 
ically, the Wilson coefficients /„ are series in ar s and en- 
code the ultraviolet physics that has been integrated out 
from QCD. The operators 0^ of dimension d„ describe 
the low-energy dynamics and are counted in powers of 
v. Heavy quarkonium annihilations are controlled by the 
imaginary part of the NRQCD Hamiltonian, i.e. the imag- 
inary part of the Wilson coefficients of the 4-fermion op- 
erators (0^" } - ^ Knxx* K' n ij/) in the NRQCD Lagrangian. 
The wave function of the quarkonium state is given by a se- 
ries of terms in which the leading one is the quark antiquark 
in a color singlet state and the first correction, suppressed 
in v, comes from quark-antiquark in an octet state with a 
gluon: 

\H) = QQQ 1 ) + \QQ s g) + ...)®\nljs). (3) 

To calculate physical quantities like spectra and decays the 
operators have to be evaluated over the wave functions and 
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the power counting of the two combines to give the order 
in v of the calculation. It is then clear that in the case in 
which the octet operators are enhanced with respect to the 
singlet, the octet part can be as relevant as the singlet one. 
The EFT contains naturally octet contributions. 

3.1 Spectrum 

The NRQCD Lagrangian is well suited for lattice evalu- 
ation. The quark propagators are the nonrelativistic ones 
and since we have integrated out the scale of the mass, the 
lattice step used in the simulation may be a factor 1 /v big- 
ger. Lattice evaluation of heavy systems like bottomonium 
become thus feasible. The latest results for the spectra 
(quenched and unquenched) are given e.g. in r il5l . The 
radial splittings are accurate up to order 0{a s v 2 ) while fine 
and hyperfine splittings are accurate only up to 0(a s ), due 
to the fact that only tree level matching coefficients have 
been used. A calculation of the NRQCD matching coef- 
ficients in the lattice regularization scheme is still missing 
and would be relevant to improve the precision of the lat- 
tice data. 

3.2 Decays 

NRQCD gives a factorization formula for heavy quarko- 
nium (H) inclusive decay widths into light hadrons (LH) [ 

El 

r(H -> LH) = ^ {U^K^Km). (4) 

n 

Similar formulas hold for the electromagnetic decays. The 
4-fermion operators are classified with respect to their ro- 
tational and spin symmetry (e.g. 0( 2S+l S j), 0( 2s+1 Pj), ...) 
and of their colour content (octet, (9 8 , and singlet, 0\, op- 
erators). Singlet operator expectation values may be eas- 
ily related to the square of the quarkonium wave functions 
(or derivatives of it) at the origin. Octet contributions re- 
main as nonperturbative matrix elements of operators over 
the quarkonium wave functions. According to the power 
counting of NRQCD, the octet contribution (h\0&C S o)\h) 
to F-wave decays is as relevant as the singlet contribution 
r 1141 . This octet contribution reabsorbs the dependence 
on the infrared cut-off of the Wilson coefficients solving 
the problem that arised in the color singlet potential model. 
Systematic improvements are possible, either by calculat- 
ing higher-order corrections in the coupling constant or by 
adding higher-order operators. If one goes on in the expan- 
sion in v, that seems to be necessary for charmonium, the 
numbers of the involved nonperturbative matrix elements 
of octet operators over quarkonium states increases in such 
a way that limits the prediction power. 

Besides this, precise theoretical predictions are also ham- 
pered by uncertainties in the NRQCD matrix elements and 
large corrections in NLO in a s . The convergence of the per- 
turbative series of the four-fermion matching coefficients is 



indeed often bad (for examples see f 1241 V A solution may 
be provided by the resummation of the large contributions 
in the perturbative series coming from bubble-chain dia- 
grams. This analysis has been successfully carried out in 
some specific cases in [ 25 1. 

3.3 Production 

As we have explained in the previous section, colour sin- 
glet production and colour singlet fragmentation underes- 
timated the data on prompt quarkonium production at Fer- 
milab by about an order of magnitude indicating that addi- 
tional fragmentation contributions were missing [ 12 1 . This 
missing contribution is precisely the gluon fragmentation 
into colour-octet 3 S i charm quark pairs. The probability to 
form a J I if/ particle from a pointlike cc pair in a colour octet 
3 S i state is given by a NRQCD nonperturbative matrix el- 
ement which is suppressed by v 4 relative to the leading sin- 
glet term but is enhanced by two powers of a s in the short 
distance part for producing color-octet quark pairs. When 
one introduces the leading colour-octet contributions, then 
the data of CDF can be reproduced. Still remains a puzzle 
the behaviour of the polarization at high p t HI 21 . 

4 Potential Non Relativistic QCD 

In NRQCD the dominant role of the potential as well as 
the quantum mechanical nature of the problem are not yet 
maximally exploited. A higher degree of simplification 
may be achieved building another effective theory for the 
low energy region of the non-relativistic bound-state, i.e. an 
EFT where only the ultrasoft degrees of freedom remain 
dynamical, while the rest is integrated out. We integrate 
out the scale of the momentum transfer ~ mv which is sup- 
posed to be the next relevant scale. Then, two different 
situations may exist. In the first one, mv » Aqcd and the 
matching from NRQCD to pNRQCD may be performed in 
perturbation theory, expanding in a s . In the second situa- 
tion, mv < Aqcd, the matching has to be nonperturbative, 
i.e. no expansion in a s is allowed. Recalling that r _1 ~ mv, 
these two situations correspond to systems with inverse 
typical radius smaller or bigger than Aqcd, or systems re- 
spectively dominated by the short range or long range (with 
respect to the confinement radius) physics. 

4.1 mv » Aqcd 

The effective degrees of freedom are: QQ states (that can 
be decomposed into a singlet and an octet wave function 
under color transformations) with energy of order of the 
next relevant scale, Aqcd, mv 2 and momentum p of order 
mv, plus ultrasoft gluons A jU (R, t) with energy and momen- 
tum of order Aqcd,'mv 2 - All the gluon fields are multi- 
pole expanded (i.e. expanded in r). The Lagrangian is 
then an expansion in the small quantities p/m, 1/rm and 
<9(A Q cd, mv 2 ) x r. 
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The pNRQCD Lagrangian is given at the next to leading 
order in the multipole expansion by I" I17I : 



^pNRQCD 



5o _iL_y j(r) _y^ 

m m" 



n=\ 



+o 



iD 



+gV A (r)Tr {o f r ■ E S + S + r ■ E o} (5) 
+^-^Tr {O'r • E O + O 1 Or ■ e} - ^F a ^ va . 

At the leading order in the multipole expansion, the sin- 
glet sector of the Lagrangian gives rise to equations of mo- 
tion of the Schrodinger type. The two last lines of con- 
tain (apart from the Yang-Mills Lagrangian) retardation (or 
non-potential) effects that start at the NLO in the multipole 
expansion. At this order the non-potential effects come 
from the singlet-octet and octet-octet interactions mediated 
by an ultrasoft chromoelectric field. 

Recalling that r ~ l/mv and that the operators count like 
the next relevant scale, 0(mv 2 , Aqcd), to the power of the 
dimension, it follows that each term in the pNRQCD La- 
grangian has a definite power counting. From the power 
counting e.g., it follows that the interaction of quarks with 
ultrasoft gluons is suppressed in the Lagrangian by v ( by 
gv if mv 2 » Aqcd) with respect to the LO. 

The singlet and octet potentials are well defined objects 
to be calculated in the perturbative matching. In this way 
a determination of the singlet QQ potential at three loops 
leading log has been obtained in [ 19 1 and consequently 
also a determination of ay which shows how this quantity 
starts to depend on the infrared behaviour of the theory at 
three loops. 

4.1.1 Spectrum 

Given the Lagrangian in Q it is possible to calculate the 



(n)\ 



O 



quarkonium energy levels at order maf [ [19] 1201 . At this 
order the energy E„ of the level n receives both from the 
average value of the potentials and from the singlet-octet 
interaction (retardation effect) a contribution that read 

En = {n\VM\n)-ig 2 dt {n^'^Mn) (E(t)E(0)>Qi)(6) 
Jo 

being H$ abd H„ the singlet and octet Hamiltonian respec- 
tively. The nonlocal electric correlator is a nonperturba- 
tive object dominated by the scale Aqcd, (E(?) E(0))(/u)) ~ 
exp{-AQCDf}- Thus the integral in (|6j is the convolution 
of two exponentials with exponent of order mv 2 (for the 
energy) or Aqcd (for the correlator). Depending on the 
relative relation of the two scales three different situations 
take place: 



• if mv 2 » Aqcd then the correlator reduces to 
the local gluon correlator and the second contribu- 
tion in (|6) corresponds to the previously mentioned 
Voloshin-Leutwyler sum rule contribution HI II . 

• if mv 2 <k Aqcd then the energy exponential can 
be expanded and the second contribution in cor- 
responds to a short range nonperturbative potential 
corrections 1 T17I . 

• if mv 2 ~ Aqcd then neither exponentials can be ex- 
panded and the nonlocal condensate has to be used 
in the energy level calculation F II 71 . 

Both the potential model and the sum rule results are con- 
tained in pNRQCD as different kinematical limits. What 
appeared as a puzzle and a problem at the origin is now 
understood as a consequence of the richness of the quarko- 
nium dynamics and is appropriately accounted for by the 
EFT. 

The calculation of the quarkonium energy levels at higher 
orders in perturbation theory is relevant to extract the 
masses of the heavy quarks from the T{IS), J/tfr and ttbar 
production cross section cf. \ 1261 1201 . The perturbative de- 
termination of the levels, have been used in [ E| [28 1 for 
the calculation of the energy levels of some lowest reso- 
nances of bottomonium, charmonium and B c , after having 
removed the renormalon (between the pole mass and the 
singlet static potential) and under the assumption (or to test 
the assumption) that the nonperturbative corrections, in the 
form of nonlocal condensates or short range nonperturba- 
tive potentials, are in total small. 

4.2 Aqcd ~ mv 

In this case the (nonperturbative) matching to pNRQCD 
has to be done in one single step [ 1211 . Under the cir- 
cumstances that other degrees of freedom (like those as- 
sociated with heavy-light meson pair threshold production 
and heavy hybrids) develop a mass gap of order Aqcd the 
quarkonium singlet field S remains as the only dynami- 
cal degree of freedom in the pNRQCD Lagrangian, which 
reads [QUEUES 



^pNRQCD 



(7) 



In this regime we recover the quark potential singlet model 
from pNRQCD [ E2 E3- The final result for the po- 
tentials (static and relativistic corrections) appears factor- 
ized in a part containing the high energy dynamics (and 
calculable in perturbation theory) which is inherited from 
the NRQCD matching coefficients, and a part containing 
the low energy dynamics given in terms of Wilson loops 
and chromo-electric and chromo-magnetic insertions in the 
Wilson loop f 1211 . Such low energy contributions can be 
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calculated on the lattice [ 30 1 or evaluated in QCD vacuum 
models [0 I31I . The expression obtained for the poten- 
tial is the QCD expression, in particular all the perturbative 
contributions to the potential at the hard scale are correctly 
taken into account. This solves the problem of consistency 
with perturbative one-loop calculations that was previously 
encountered in the Wilson loop approach. Moreover, fur- 
ther contributions, including a 1/m nonperturbative poten- 
tial, appear with respect to the Wilson loop original results 

IE- 

4.2.1 Decays 

The inclusive quarkonium decay width achieve in pN- 
RQCD a factorization with respect to the wave function 
(or its derivatives) calculated in zero which is suggestive 
of the early potential models results: (cf. eq. Q) 

T(H -> LH) = F(a s , Aqcd) • l<A(0)| 2 . (8) 

Similar expressions hold for the electromagnetic decays. 
However, the coefficient F depends here both on ff s and 
Aqcd- In particular all NRQCD matrix elements, including 
the octet ones, can be expressed through pNRQCD as prod- 
ucts of universal nonperturbative factors by the squares of 
the quarkonium wave functions (or derivatives of it) at the 
origin. The nonperturbative factors are typically integral 
of nonlocal electric or magnetic correlators and thus de- 
pending on the glue but not on the quarkonium state H22I . 
The presence of this nonperturbative correlators is indi- 
cated by the Aqcd dependence of F in Typically F 
contains both the NRQCD matching coefficients / at the 
hard scale m and the nonperturbative correlators at the low 
energy scale Aqcd- The nonperturbative correlators, being 
state independent, are in a smaller number than the non- 
perturbative NRQCD matrix elements and thus the predic- 
tive power is greatly increased in going from NRQCD to 
pNRQCD. In [ 1221 the inclusive decay widths into light 
hadrons, photons and lepton pairs of all S -wave and P- 
wave states (under threshold) have been calculated up to 
0(mv 3 x (AQ CD /m 2 , E/m)) and 0(mv 5 ). A large reduc- 
tion in the number of unknown nonperturbative parameters 
is achieved and, therefore, after having fixed the nonper- 
turbative parameters on charmonium decays, new model- 
independent QCD predictions are given for the bottomo- 
nium decay widths [ 22 1. 

Once the methodology to compute the potentials (real and 
imaginary contributions) and from these the inclusive de- 
cays, within a 1/ra expansion in the matching has been de- 
veloped, the next question to be addressed is to which ex- 
tent one can compute the full potential within a 1/m expan- 
sion in the case Aqcd s> mv 2 . It has been shown [ 23 1 that 
new non-analytic terms arise due to the three-momentum 
scale sjm Aqcd- These terms can be incorporated into lo- 
cal potentials (£ 3 (r) and derivatives of it) and scale as half- 
integer powers of 1/m. Moreover, it is possible to factor- 
ize these effects in a model independent way and compute 



them within a systematic expansion in some small param- 
eters [ 23 1. 

4.2.2 Production 

Since the power counting of pNRQCD may be different 
from the power counting of NRQCD, we expect that we 
may eventually explain in this way some of the difficulties 
that NRQCD is facing in explaining the polarization of the 
prompt J I iff data. In particular, if the magnetic field turns 
out to be not suppressed with respect to the electric field 
operator in the power counting, then the spin flip term is 
enhanced and the polarization may be diluted explaining 
the behaviour of the data with high p, r 1221 |2"T1 1291. 

4.3 Renormalization group improvement and 
Poincare invariance 

The effective field theory can be used for a very efficient 
resummation of the large logs (typically logs of the ratio 
of energy and momentum scales) once a renormalization 
groups analysis of the EFT has been performed. Such pro- 
gram has been successfully performed in pNRQCD f 1321 . 

Since the EFTs are constructed to be equivalent to QCD, 
Poincare invariance has still to hold. In [ 33 1 the constraints 
induced by the algebra of the Poincare generators on non- 
relativistic effective field theories have been discussed. It 
has been shown that Poincare invariance imposes well de- 
fined relations among the EFT matching coefficients. The 
relations have been given both for NRQCD and for pN- 
RQCD. 



5 Conclusions 

The progress in our understanding of non-relativistic effec- 
tive field theories makes it possible to move beyond ad hoc 
phenomenological models and have a unified description 
of the different heavy-quarkonium observables, so that the 
same quantities determined from a set of data may be used 
in order to describe other sets. The old puzzles and prob- 
lems have been clarified and have been understood inside 
the EFT formulation. Predictions based on non-relativistic 
EFTs are conceptually solid, and systematically improv- 
able. EFTs put quarkonium on the solid ground of QCD: 
quarkonium becomes a privileged window for precision 
measurements, new physics, confinement mechanism in- 
vestigations [[Q. 
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